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Abstract 

In this paper we consider random Hamiltonians defined on long-range percolation 
graphs over II 1 . The Hamiltonian consists of a randomly weighted Laplacian plus a 
random potential. We prove uniform existence of the integrated density of states and 
express the IDS using a Pastur-Shubin trace formula. 



on 

1 Introduction 

In the study of solution properties and spectral features of random ergodic and periodic 
operators the integrated density of states (IDS), also known as spectral distribution function, 
■ plays and important role. While it is by its very nature a much simpler object than the 

original one, i.e. the operator family, it exhibits a number of features which turn out to be 
crucial for the understanding of the spectrum and the corresponding eigensolutions. Let us 



(S) 1 spell out some of these features explicitly: 

o 



■ (A) the low energy asymptotics of the IDS, 



(B) its local and global continuity properties, and 

(C) its approximability by finite volume analogues. 

X 

In fact, (C) is of relevance for all other questions about the IDS since all methods used 
to answer /understand them rely on one stage or another on finite volume approximations. 
Of course, depending on the question which is being considered, the type and quality of the 
approximation (C) will vary. 

Let us highlight the intimate relationship between the continuity (B) and approximability 
(C) properties on an elementary level, to provide a motivation for the discussion which follows. 
On the one hand, if a sequence of probability measures converges weakly to a measure with 
no atoms, the corresponding distribution functions converge already uniformly. On the other 
hand, uniform convergence of continuous distribution functions obviously implies that the 
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limiting measure has no atoms. Now, there are important classes of random operators where 
it is known that the IDS is not continuous (and other, where continuity is still an open 
question). In this setting the question arises whether the approximability (C) in the L°°- 
topology persists or not. More precisely, it is natural to ask, what are reasonable criteria, 
and which is the mechanism, which ensures that the finite volume approximating distribution 
functions converge uniformly to the IDS. 

Let us review the history of works addressing the above explained approximability ques- 
tion. In the simplest setting of finite-hopping range ergodic operators on £ 2 (Z, d ) the continuity 
of the IDS was established already in |DS84j . In fact, this result can be sharpened under very 
mild conditions to a quantitative form of continuity, namely log-Holder continuity, c.f. [CS83aJ. 
In the corresponding continuum setting, more specifically for ergodic Schrodinger operators 
on L 2 (M. d ), the complexity of this question is highly dependent on the space dimension d. 
For one dimension, log-Holder continuity of the IDS of ergodic operators was established in 
|CS83bj . This result uses a technique very specific to one space dimension, namely trans- 
fer matrices. The higher dimensional case was tackled only very recently in [BKllj . where 
log-Holder continuity of the IDS was shown for Schrodinger operators on L 2 (M. d ) for d < 3. 
The case d > 4 is still open. This is due to the fact that the present state of knowledge 
about quantitative unique continuation properties of eigensolutions of Schrodinger operators 
based on Carleman estimates, cf. e.g. |BK05j . is good enough for the IDS-continuity proof 
in dimensions three and less, but not above this threshold. (Let us also note that this result 
is not stable under the addition of magnetic field, the primary counterexample being the 
Landau Hamiltonian, whose IDS is a step function. Note that the eigenfunctions there have 
a Gaussian decay, in contrast to the lower bounds on the decay established in [BK05] .) 

Once one moves away from the Euclidean setting, continuity is no longer a prevail- 
ing feature of the IDS. There are two important prototypical examples of discrete Laplace 
and Schrodinger operators whose IDS exhibit jumps: Hamiltonians on percolation clusters 
jCCF+86] . [Ves05j and quasi-crystal graphs [LS031 ILS061 IKLS03] . In the former model the 
discontinuities form even a dense set in the spectrum. Nevertheless, it has been established 
that the sequence of normalized finite volume eigenvalue counting functions converges at every 
energy and even uniformly along the energy axis to the IDS, cf. |Ves051 [LS061 ILMV081 ILV09] . 

The question whether the approximability (C) holds at any given energy is also of in- 
terest in a purely geometric setting, where no randomness in involved. More specifically, 
the papers [Liic94j IMY02} IMSY031 lDLM+03] study the approximation of Betti numbers by 
their finite volume analogues. The usual weak convergence of measures is not sufficient to 
yield this statement. The results of [LV091 ILSVIOI |PS12| actually apply under very mild and 
natural geometric assumptions (namely amenability) and thus lift the pointwise everywhere 
convergence results of |Luc94l IMY021 IMSY031 lDLM+031 to uniform convergence. 

However, all papers mentioned so far concern local operators (in the continuum setting) 
respectively finite-hopping range operators (in the discrete setting). Let us spell out the last 
property explicitly: A bounded operator A on £ 2 (G) for some graph G is called of finite 
hopping range R S N if for any tp £ £ 2 (G) and x € G with distance larger than R to suppt/? 

{A<p)(x) = 

The first result for an ensemble of operators on £ 2 {G) beyond this restriction was achieved 
in [Schl2j . It studies Laplace operators on long-range percolation graphs. These operators 
are not of finite hopping range, albeit non-zero "matrix elements far off the diagonal" appear 
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with small probability. Due to the methods applied in |Schl2j , the randomness present in the 
model had to be of finite local complexity. 

In the present paper we continue this line of research. Here we are able to treat models 
with real-valued entries, possibly continuously distributed, and with long-range interactions. 
As mentioned above, the proof of |Schl2] does not apply in this setting, but a combination 
of ideas from this paper together with methods from |LV09| allows one to prove uniform 
approximability in this more general situation. Also, we will treat here randomized versions 
of adjacency as well as Laplace operators. The results of |Schl2| apply to the second type of 
ensembles only. 

Question (A) has been addressed for certain long-range percolation models before in 
|AV09j . extending previous results for the usual quantum percolation model on Z d , cf. for 
example [KM06, MS07 . We will quote for completeness sake the result of |AV09j below, once 
we have the necessary notation at disposal. 

Let us stress an important feature of long-range percolation Hamiltonians. They provide 
a simple model interpolating between discrete random Schrodinger operators and random 
matrices. This is of interest, since the two last mentioned classes of operators have quite 
different spectral features. Thus one is led to ask in which aspects and regimes, long-range 
percolation Hamiltonians share features with one or the other of these classes. To explain the 
structural difference between discrete random Schrodinger operators and random matrices 
let us restrict ourselves for the moment to operators on finite segments of Z. The arising 
Schrodinger operators are a special type of finite Jacobi matrices, and are in particular tri- 
diagonal. Random matrices have a full array of random entries. The non-zero entries of a 
Laplacian of a typical long-range percolation graph are concentrated mostly near the main 
diagonal, however no diagonal will consist entirely of zeros. This is the mentioned interpolat- 
ing property. The relation between long-range percolation models and random matrices was 
for instance studied in |Aya09a Aya09b| . 

Let us describe the content of the paper in detail. In the next section we state the main 
result in a concise form, discuss extensions to groups and graphs beyond Z d and the relation 
to results concerning the low energy asymptotics of the IDS (A) established for long-range 
percolation graphs in [AV09J. 

In Section [3] we present the long-range percolation model and in particular the underlying 
probability space. Furthermore we define the randomly weighted Hamiltonian H u on the long- 
range percolation graph as a selfadjoint and metrically transitive random operator. Depending 
on the choice of the corresponding parameters, the operator in question is either the adjacency 
operator of a long-range percolation graph, a Laplacian, or a Schrodinger operator (i.e. a 
Laplacian plus a random potential), each one with random weights on the edges. For a 
realization oj the restriction of this operator to a finite box A n C Z d gives rise to the eigenvalue- 
counting function F£ . This function encodes the distribution of the spectrum of the restricted 
operator. In the sequel we pursue the question whether and in which sense the limit of 
the sequence of functions (F£) exists. In Section HI weak convergence of this sequence is 
established for almost all realizations using a result of Figotin. In Section [5] this statement 
is upgraded to uniform convergence along the energy axis. This is done by proving that the 
functions F% do not only approximate the limit function at its continuity points, but also 
give an efficient estimate of the size of the jumps at a point of discontinuity. The key tool to 
control the size of the jumps is Theorem 1 5. 5 \ formulating the main technical contribution of 
the present paper. 
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Let us stress that we give a detailed account of all the main steps of the proof thus making 
it accessible to non-specialists. 



2 Main result 

Here we formulate the main result of the paper using a minimum of notation needed for this 
purpose. More detailed definitions of the framework can be found in Section [3j 

Denote by E := {{x,y} C Z d | 1,1/ £ Z d } the set of all edges (or loops) on 7L d and 
by a e ,b e ,e € E a collection of independent real- valued random variables on a probability 
space (fl,A, P). The variance of the <i e ,e 6 E is uniformly bounded and any two &{ x ,y} &nd 
a {x+k,y+k} are identically distributed for x,y,k € Z d . Let p € £ 1 ('Ii d ) with < p(x) = p(—x) < 
1 and 6{ XjJ/ } be Bernoulli distributed with parameter p(x — y). The kernel, resp. the matrix 
entries of the random Hamiltonian H(u) are given by: 

H m M .= HSf( U ) .= if (21) 

[aa {x} {u)b {x} (u) - f3 Y.z+x a {x, z }{u)b {x ^{u) ifx = y, 

where a, ft G [0, 1] are fixed numbers. Selfadjointness properties of H(lj) are discussed in 
Section [3l Depending on the values a and (3 we obtain several interesting subclasses: ran- 
domly weighted Laplacians or adjacency operators, with or without random potentials (on 
the diagonal), cf. Remark 13.41 

For neN let A n := ([-n,n] n Z) d , H n (uS) be the restriction of H(u) to A ra , and F%(\) 
the number of eigenvalues of H n (uj) not exceeding A, counting eigenvalues according to their 
multiplicity. Set F : R -> R, 

F(X) := E{(E H ((-oo, \})5 , S )}, (2.2) 
where E H ^((— 00, A]) is the spectral projector of H(uS). 

Theorem 2.1. Lei F£,F be the distribution functions given above. Then there exists a set 
Q C U of full measure such that for all u G O we /iaue 



lim sup 

n— >oo 



^-AX)-F(X) 



0, 



I A, 

This gives a representation of the IDS by a Pastur-Shubin trace formula. 



Extension to more general geometries 

In the presentation of our results we have not striven for the maximal possible generality, 
but rather tried to present proofs as explicitly as possible in their most accessible form. This 
concerns in particular the restriction to operators defined on £ 2 (7j d ). In fact, the presented 
results carry over to much more general graphs T (more precisely to operators on £ 2 (F)). 

We leave the explicit calculations for another occasion, cf. jSchj . and state only the scope 
of the general results: 

Let G be a finitely generated discrete amenable group, and T a graph on which G acts freely 
and cocompactly by translations. Completely analogously as in the Z rf setting one can define 
long-range percolation on such graphs, as well as random operators introduced in Definition 
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(|2.ip . respectively (|3.5p . For such models the results which we use from [PF92j can be proven 
analogously. Consequently, Lemma 13.61 and Theorem 14.21 have their generalized counterparts 
in this setting. Relying on the ideas of |LV09j and |Schl2j one can see that it is possible 
to extend the results of Section as well. Here one needs to apply the pointwise ergodic 
theorem of Lindenstrauss [LinOl] instead of Theorem 15.11 to finally obtain that Theorem 12.11 
holds analogously. 

Let us note that the method we use here is in the sense efficient, that it does not need 
any condition beyond amenability, i.e. the existence of a F0lner sequence, on the discrete 
group G. In comparison to this, the method of |Schl2j relies on the following additional tiling 
condition: 

It is assumed that there exists a F0lner sequence (Q n ) n such that for each n € N there is 
a set T n = T~ l C G with the property that G is the disjoint union of the sets {Q n t \ t G T n }. 
This assumption is satisfied for many amenable groups, however it is not clear whether it 
holds for all of them. 

Low energy asymptotics 

In [AV09j the authors study low energy asymptotics of the IDS for the long-range percola- 
tion model. The operator under consideration is the (non- weighted) Laplacian A w , i.e. the 
operator acting on any ip: T, d — > C with finite support by 

(A u <p)(x) = ]T {<p(x) - <p(y)) 

This equals the operator in (|2.ip in the case a = 0, ft = 1 and where the random variables a e 
are constant 1, see also (|3.6|) 

An in physical communities common way to introduce the long-range percolation model 
is the following: For each pair of vertices x, y € Z d let J XjV be a real number such that 

• Jx,y — Jy,x 

• J x +z, y +z = Jx,y for all z G Z d , 

• J := J x := J2y£Z d Jx,v ls finite and independent of x € 

We fix ft > and declare an edge {x, y} to be an element of with probability 1 — e~@ Jx ' y . 
This gives the random graph = (Z, d ,E UJ ). Notice that the probability that certain edge 
is an element of E u is increasing in f3. Thus, the subcritical phase, in which all clusters 
are almost surely finite corresponds to small values of the parameter j3 and the supercritical 
phase in which there exists almost surely an infinite cluster corresponds to large values of the 
parameter ft. Just like in the case of the nearest neighbor percolation model these two phases 
are separated by a single value of the parameter ft. The authors of |AV09| define the IDS as 
in (|2.2p and prove that for every subcritical ft there are constants c(ft),d(ft) > such that for 
E > small enough 

exp (-cO^- 1 / 2 ) < F(E) - F(0) < exp (-d{ft)E~ 1 / 2 ^ . 

Actually the results of [AV09] apply to operators on quasi-transitive graphs Y. The present 
result is complementary to these observations, as we show that the finite volume approximants 
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do actually converge to this limit function given by the Pastur-Shubin formula. Furthermore 
the combination of both results shows that even the approximating functions will exhibit 
exponential behavior for low energies. 

Note that in this paper we introduce the long-range percolation model in another, but 
equivalent (see |Schl2j ) way, via a certain function p G see (|3.ip . More background 

on the models considered in [AV09| can be found in the review paper |AV08j . 

3 Setting and first results 

Let r be the Z rf lattice and denote by d : Z d x 7L d — > No the graph distance in the lattice 
or equivalently the ^-distance in Z rf . With this metric we define the i?-boundary of a set 

A C 7L d by 

d R A := {x G A | d(x, y) < R for some y G Z d \ A}. 

Furthermore we let E := {{x, y} C 7L d \ x, y G Z rf } be the set of all subsets of Z d containing 
either one or two elements. The set E can be interpreted as the edge set of the complete 
undirected graph over Z rf , containing loops at each vertex. 

The probability space (0, A, P) is given in the following way. The sample space is $7 = 
f[ egE (M x {0,1}) and we denote the elements in by uj = (uj' e ,uj^) eG E- The appropriate 
<7-algebra is A = &) eeE (B(R) ® ^({A I}))- I n order to define a measure on this space we fix 
some p G £ 1 (7 J d ) with 

< p(x) = p(-x) < 1 (x£Z d ) (3.1) 
and for each z G Z rf some probability measure \i z on R such that there is v € R with 

/ x 2 dfi z (x) <v 2 (z G Z d ). (3.2) 
Jr 

We set P := <S>{ Xjy y eE (^ 

-y <8) v x-y) where for each z € 7L d the measure v z is Bernoulli with 

parameter piz). 

Remark 3.1. The cr-algebra A is generated by the cylinder sets which are given the following 
way 

Z = {Z(A ei ,B ei ,...,A ek ,B ek ) \ ken, ei eE,A et eB(R),B £j £?({0,1}) for i = 1, . .. , k} 
where 

Z(A ei ,B ei ,...,A ek ,B ek ) = {uj G n | J e . eA ei ,oo". G B ei for t = 1, . . . , k} . 

Now for each uj = (uj' e ,uj") e< =E and e G E we set a e (o;) := and 6 e (w) := a;". This 
procedure gives independent random variables a e ,b e , e G E satisfying P(a e G B) = fi e (B) 
as well as P(6 e = 1) = u e ({l}) = p(x — y) for arbitrary e = {x,y} G E and G £>(R). 
Furthermore, by (j3.2|) we have for each e G E 

E(|a e |) <v 2 + l. 

These random variables induce for each oj G a graph T w = (Z^,^) with weighted edges. 
Here Z d is the vertex set and E w is the subset of E, where an edge e € £ is an element of E w 
if and only if b e (oj) = 1. In this case one can think of a e (uj) a the weight of the edge e. For a 
subset A C 7L d and an element x G Z rf we write x ~ A if there exists y G A with {x, y} G E u . 

The following Lemma shows that T w is almost surely locally finite, i.e. each vertex is 
incident to only finitely many edges in T u . 
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Lemma 3.2. The graph is locally finite for almost all wed. 

Proof. Fix an element x G 7L d and consider the events A y := {bi Xiy \ = 1}, y G Z d . Then 
clearly 

£ P(^) = ^ p(x - y) < oo, 

as p G Hence, the Borel-Cantelli Lemma gives a set Cl x of full measure such that each 

oj G Oj; is contained in only finitely many A y , y G As Z rf is countable £7 := f] x ^i d ^-x is a 
set of full measure as well. Furthermore is locally finite for all oj E Cl. ■ 

Given 7 G Z rf , let us define translations T 7 : Q, — > by 

T 7 (w) = T 7 ((Wg,Wg) eg£ ;) = (Wg +7 ,Wg +7 ) eg£ ; 

where for e = {5, /i} € £ we mean by e + 7 the element {g + 7, /i + 7} G E. For 7 G Z rf and 
i? G .4 we denote the image and the preimage of B under T 7 by 

T y (B) = {T 7 (w) G fi I w G B} and T^B) = {w G fi | T 7 (w) G 5}. 

Note that for B G i we have T" 1 ^) = T y -i(B). We further define T to be the mapping 
7 1— > T 7 which maps each element of into the space of automorphisms on (fl,A,¥). Note 
that by definition T is ergodic if and only if for any B G A with T^(B) = B for all 7 G Z rf 
one has P(.B) G {0, 1}. The following result is basic, but we do not know an explicit reference 
in the literature, so we include a proof for completeness sake. 

Lemma 3.3. T is a measure preserving, ergodic left-action on (Q,A,¥). 

Proof. For an edge e = {g, h} G E, vertices x, y G Z rf and wsfiwe have Tq(oj) = uj and 

T x+y (uj) = (uj' e+x+y ,u}" +x+y ) e <z E = T x (T y (u)) 

which shows that T is a left action of Z d on fL 

By definition of P and the random variables a e and b e we have P(a e G B) = P(a e + 7 G B) 
as well as P(6 e = 1) = P(6 e+7 = 1) for any e G E, 7 G 7L A and B G B(R). Furthermore, as T 7 
is a translation, P(Z) = P(T 7 (Z)) holds obviously for any 7 £ Z d and any cylinder set Z G Z, 
which implies the same property for any set B G A, c.f. Remark 13.11 

To prove ergodicity let B E A with B = T^{B) for all 7 G Z d and F(B) > be given. 
We need to show that this implies ¥(B) = 1. In the following we apply the approximation 
lemma for measures, which belongs to the entourage of Caratheodory's extension theorem, 
cf. e.g. Theorem 1.65 in |Kle08j . Let e > 0. As B G A = a{Z) and as Z is semiring we can 
find cylinder sets Z\ , . . . , Z n E Z such that 

n 

F(BAZ) < e where Z := (J Z k , 

k=l 

which gives 

F{B) 2 - 2F(B)e < F(Z) 2 < ¥(B) 2 + 2F(B)e + e 2 . (3.3) 
Furthermore we have for any 7 G Z d 

F(Z n T 7 Z) <P((BU(Z\ B)) n T 7 Z) 
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<P(Bn T 7 Z) + P ((Z \ B) n r 7 z) 
<P(Bn (T 7 B U (T 7 Z \ T 7 S))) + e 

<p(Bn r 7 s) + p(en (r 7 z \ r 7 B)) + e 

<P(Bn T 7 S) + 2e 
By symmetry we get for all 7 G Z d 

P(5 n T 7 5) - 2e < P(Z n T 7 Z) < P(B n T 7 S) + 2e 
and the T-invariance of B implies 

¥(B) -2e< P(Z n T 7 Z) < P(J3) + 2e. 



(3.4) 



As Z is a finite union of cylinder sets, it does only depend on finitely many edges. Hence 
there exists an element h G Z d such that Z and T^Z are independent, which gives 

F(Z n T h Z) = F(Z)F(T h Z) = F(Z) 2 



since T is measure preserving. This gives together with (|3.3p and f|3.4|) 

P(S) - 2P(B)e - e 2 - 2e < F(B) 2 < F(B) + 2¥(B)e + 2e 
and dividing by P(B) > leads to 

1 _ 2e - < F(B) < 1 + 2e + 



As these inequalities hold for arbitrary P(£?) > e > we get P(2?) = 1. ■ 

Denote by £ 2 (Z d ) all square summable, complex-valued functions on 7L d and by C c (Z rf ) the 
subset of £ 2 (Ij d ) consisting of all finitely supported functions. Let a,/3 G [0, 1] be some fixed 
numbers. Using the random variables a e ,b e , e G E we define for each uj G VL as in Lemma 
the random operator H(u) := H a ^(uS) : C c (Z d ) ->■ £ 2 (Z d ) point-wise by 

fl* v (w) := := (^J W ^ (U) lf X + V > 

\aa {x} {uj)b {x} {uj) -^E^ a {x,z}(w)^, 2 }H if x = y 



and for G C c (Z d ) we set 

(^(w)^)(x) := (i^(u,)^)(x) == E H*,v(foMv)- 



(3.5) 

It is easy to see that 

(H(u)(p)(x) = (<P(v) - P<f( x )) a {x, y }(u) + aip(x)a {x} (uj). (3.6) 

Using this we obtain for each ip G C C (Z 2 ) and uj G such that is locally finite that 
H(u)tp G ^ 1 (Z d ) C ^ 2 (Z d ). To see this we set ^4 := suppy?, m := max^g^ |y(^)| and N y (uj) := 
{x G Z d I {x, y} G -E^} to estimate 
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y€Ax£N v (oj) y&A 



\Hy,y(u)\ + V \H Xty (uj] 



V 



xeN y (u) 




< oo. 



Note that here we used that N y (u}) is finite, as the underlying graph is locally finite. 
The mapping H : f2 — > L(^ 2 (Z d )), w i— > H(u) is a random operator with domain C c (Z d ), 
c.f. [PF92] , Here L(£ 2 (Z d )) is the space of the linear operators which are densely defined in 
^ 2 (Z d ). 

Remark 3.4. The operator H{to) depends on the choice of a, (3 € [0, 1] and is defined on the 
finitely supported functions in £ 2 (7j d ). In Lemma 13.61 we will define the self-adjoint extension 
H(u) of this operator. Depending on a and (3 we have in particular the following special 
cases for H(uj): 

• if a = and /3 = 1, then H(uj) is the randomly weighted Laplacian on the graph r w , 

• if a = f3 = 1, then H(u>) is the randomly weighted Laplacian on the graph T w plus a 
random diagonal, 

• if a = 1 and /? = 0, then H{uj) is the randomly weighted adjacency operator of T w with 
a random diagonal, 

• if a = /3 = 0, then H{uj) is the randomly weighted adjacency operator of T w with zeros 
on the diagonal. 

The diagonal elements which appear if a > can be interpreted, either as random weights 
on the loops or as a random potential. For values a, (3 € (0, 1) the operator can be seen as 
an interpolation between, the adjacency operator and the Laplacian respectively Schrodinger 
operator of the graph T^. 

We will use the same symbol T 7 for a mapping T 7 : L{l 2 {Wj d )) — > L(£ 2 (Z, d )) defined by 

T 1 ((A Xj y) x y€li d) := {A x+1 ^ y+1 ) x ^ y&Ij d, 

for arbitrary A = {A x>y ) x>y&d G L(£ 2 (Z d )). We set XJ 1 : £ 2 (Z d ) -> £ 2 (Z d ) 

U-y{{<-p{x)) x& d) := {(p(x + i)) x& i 

where <p = (f(x)) x£Z d is arbitrary. Then obviously T^(A) = U^AU^ 1 . 

For each x, y, 7 € Z d with x ^ y and = (Wg,w") ee f; we set s := {x, y} and have 

H x>y (T 7 (u))) = a s (T 7 (o;))6 s (T 7 (w)) 

= a s ((a;g +7 , Wg +7 ) ee £;)6 s ((a;g +7 , Wg +7 ) ee _B) 
= £^ +7 • cj" +7 

~~ W {x+7^+7} ' UJ {x+j,y+-/} 

= a {x+7,y+7}( a; ) ' ^{x+7,j/+7}( a; ) = Hx+-y,y+-y ( w ) 
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Furthermore we have for the diagonal elements 

(T 7 (a;))6 w (r 7 M)-/3^ 

a {x,z}( T l(^))b{x, Z }( T "/(^)) 

z^x 

= aa {x+ 1 }( UJ ) b {x+ 1 }( UJ ) - P'^2 a {x+j,z+ 1 }(, UJ )b{x+ 1 ,z+j}(, UJ ) = H x+1 ^ x+1 (uj). 

z^x 

Therefore we have 

H{T 1 {u))=T 1 (H{u:)). (3.7) 

Definition 3.5. Let A be a random operator mapping each element of the probability space 
(SI, A, P) to an linear operator on the Hilbert space %. Then A is called metrically transitive, 
if there exists a group T of measure preserving automorphisms of (SI, A, P), a group of unitary 
operators IA := {Ut \ T G T} on % and a homomorphism from T to U such that 

B G A such that TB = B for all T G T =4> F(B) G {0, 1} (3.8) 

and one has for all uj G SI and all T G T the relation 

A(Tu) = U T A{u)U^ x . (3.9) 

The next aim is to prove that H is essentially selfadjoint and that H and its selfadjoint 
extension are metrically transitive. 

Lemma 3.6. Let (SI, A, P) and the random operator H be given as above. Then 

(a) there exists a set SI' of full measure such that for each uj G SI' the operator H(uj) is 
essentially self-adjoint. We denote the closure of H(uj) by and its domain by D(uj). 

(b) the random operators H and 

H : n -> L(£ 2 (Z d )) given by H(uj 



I Id else 



are metrically transitive. Here Id is the identity operator in L(£ 2 (7j d )). 



To prove the Lemma we will make use of the following theorem due to Figotin Fig87 
see also [PF92| . 

Theorem 3.7 ( |Fig87| , [PF92j ). Let (Sl,A,F) be a probability space and A a metrically 
transitive random operator with domain C c (Z d ) satisfying 

E I ( \ A o,x\) ) < °°- (3-10) 

Then the operator A(uj) is for almost all uj G SI essentially self-adjoint. 
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Proof of Lemma \3.6\ First we show that H is metrically transitive. To this end, define T and 
U as follows 

T := {T 7 | 7 G Z d }, U := {U Tl := U 7 | 7 G 

and set 9? : T —> U, ip (T 7 ) = ?7 7 , which clearly is a homomorphism. It is obvious that T 
and are groups and it is easy to prove that each f/ 7 is unitary. Furthermore we know from 
Lemma [3. 31 that the translations T 7 are measure preserving automorphisms of the probability 
space. Property (|3.8p follows from the ergodicity of T shown in Lemma 13.31 as well. From 
line USD we infer that fl3jl holds. 

In order to apply Theorem 13.71 to show that H is almost surely essentially selfadjoint, it 
remains to prove (|3.10|) for H. Therefore we consider for each u> G fi 

( £ \HoA<")\) = (\Ho,o(")\ + E \ H °A")\) 

< (a\a {0} (u)b {0} (u)\ + (^ + l)^\a {oM (u)b {0>x} (u)\) ( 3 - 11 ) 



< 



4 ( E \ a {o,x}(^)\b{o, x }(^) 

x&, d 



x^O 

2 



For each uj £ ft set N(w) := {x G 7L d \ bfQ jX y(u)) = 1}. By Lemma 13 . 2 1 there exists a set & C 
of full measure such that |JV(o;)| < 00 for all uj G fi. For w G O we have 

E Ko,z}MI&{o,z}M) = ( E KtwMl) ^ l-^MI • E KtwMI 2 

a'gZ d x€N(u>) x£N(u) 

Taking the expectation value on both sides and the application of the monotone convergence 
theorem leads to 

E (( E HoM^hoM^)) 2 ) ^ E (( E ho^l 2 W^l)) ^ v 2 E E ( 6 {o,*}M)> 

^ xeZ d ' ^ x&L d ' x&L d 

where v 2 is the upper bound for the second moments given in (|3.2p . For each x G 7L d and 
w S O we set N x (oj) := |iV(w) \ then we obtain for fixed x G Z rf 

00 

E(&{o,*}M) = E fc • F ( fe {o,*}M = = fc - 1) = P(&{o,z}(") = l)E(^x + 1). 

fc=i 

Using E(N X ) < E(\N\) = \\p\\i < 00 this implies 

2 N 



+ IIpHiJ < 00. 



This shows together with (|3.11|) the finiteness of the expression in (|3.10p for the operator H. 
Hence Theorem 13.71 gives a set fi' of full measure such that for each the operator H(uj) 

is essentially selfadjoint. This proves of part (a). 

To complete the prove of part (b) it remains to show that the operator H is metrically 
transitive. This follows by the same argument as we used to prove metrically transitivity of 
H. Note that here we use that (13.711 hold for H as well. ■ 
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The operator H denned as in Lemma 13.61 is a random operator with domain C c (Z d ), c.f. 
|PF92j . We will refer to this operator as weighted Hamiltonian on the graph T u . 
Let (A n ) be a sequence of cubes given by 

A n := ([-n, n] n T) d (n G N) (3.12) 

and for each n G N let H n (uj) be the restriction of H(u) to A n . To be precise, for A C 1, d let 
the inclusion ia ■ £ 2 (A) — > £ 2 (Z d ) and the projection pa ■ £ 2 (Z d ) — > £ 2 (A) be given by 

(iA(p)(x) = ^ 6 A and (ph_ip){y) = tjj(y) 

for all x G Z d , y G A, </? G £ 2 (A) and ^ G £ 2 (Z d ). Then we set 

ff„M := p An H(uj)i An : £ 2 (A n ) -> £ 2 (A n ) 

For each w€0 and n G N we define a function i 7 ^ : R -> R by 

i^(A) := |{ eigenvalues of H n {uj) not exceeding A}| , (3.13) 

where we count the eigenvalues with their multiplicity. Therefore F£ is the distribution 
function of a measure which we will denote by pfi*\ Note that |A n | -1 p^ is a probability 
measure. 

4 Weak convergence 

In order to prove weak convergence of the approximating distribution functions we make us 
of an abstract result by Figotin |Fig87| , see also Theorem 4.8 in |PF92] . which we now cite 
in a special case. 

Theorem 4.1 ( |Fig87| , [PF92j ). Let (£l,A, P) be a probability space and A a metrically 
transitive random operator with domain C c (Z d ) such that 



]Te(LVI)<°c 



and assume that A is almost surely self-adjoint. Then there exists a set Cl C Q of full measure 
such that for all uj G f2 and all A G {s G R | -F «s continuous in s} one has 

Hm ^ = F(A) 

n->oo |A n | 

where the limit F : R — » [0, 1] given 6j/ A i — > E{(-EU(( — oo, A])5o, <5o)} a distribution function 
of a probability measure. Note that here E^ u ^((— oo, A]) is i/ie spectral projection in the 
interval (— oo, A] o/ £/ie operator A(oj) and 5 X G £ 2 (7L d ~) denotes the element with 5 x (y) = I if 
x = y and 5 x (y) = else. 

This theorem and the previous considerations immediately give the following theorem. 
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Theorem 4.2. Let the probability space (Q,A, P) and the operator H be the weighted Hamil- 
tonian given in Section 2. Set 

F : R [0, 1], F(X) := E{(£ H ((-oo, A])5 ,^o)}, 

i/ien i/iere exists a set QCfi of full measure, such that for alluj £ O i/ie distribution functions 
i^/|A n | converge to the distribution function F point-wise at all points of continuity of F. 

Proof. By definition the operator H{uj) is self-adjoint for all oj £ Vt. Furthermore H has 
domain C c (Z d ) and H is metrically transitive by Lemma [3T6l The finiteness of ^icez d ^(l^o,x|) 
follows from 

E(l^b,x|) < aE(|a {0} |6 {0} ) +/3E^|a {0i:c} |6 {0iX . } j + E(|a {0 ,*#{o,x}) 
< 2 ^ E(|a {0iZ} |& {0>x} ) < 2( V 2 + 1)11^11! < oo 

where v 2 is the upper bound for the second moments given in (|3.2p . Hence, Theorem 14.11 
implies the claim of the theorem. ■ 

5 Control of the jumps 

The aim of this section is to control the jumps of the limit function given in Corollary 14.21 in 
order to obtain uniform convergence of the approximants. In the following we will make use 
of Birkhoff's ergodic theorem in the d-dimensional case, see |Kel98] . 

Theorem 5.1. Let Z, d act from the left on a probability space (Q,A,F) by an ergodic and 
measure preserving transformation T an let (A n ) be the sequence of cubes given as in (|3,12p . 
Then for any f £ L 1 (P) 



i^oo \ A n \ J 



geA 

holds almost surely. 

Lemma 5.2. Let (Q, A, P) be the probability space, H the randomly weighted Hamiltonian 
and (A n ) be the sequence of cubes given as in (|3.12p . Then there exists a set f2 C of full 
measure such that for all oj G Q and all A £l we have 

v ^(XAn-EffMCM)) w r /TP /r UU x \-> 
LAJ = E {( E H(u)({^})So,5o)}. 

Proof. Let oj € ft be fixed. By definition of the trace we have 

Tr( XAn ^ M ({A})) = Yl (xA n E H{uj) ({X})5 x ,5 x ) = £ (E H(uj) ({X})5 x , S x ) . (5.1) 

x£Z d x£A n 

Let Fiig(H(uj), A) denote the eigenspace of H{oj) corresponding to the value A, which could 
possibly be empty if A is not an eigenvalue. Given 7 € Z d , we have ip £ Eig(i?(u;), A) if and 
only if T 7 (ip) G Eig(T y (H (u)) , A). 
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Using this we prove 

(E Tz{H{ul)) ({\})5o,5o) = (E H[oj) ({X})6 z ,S z ) . (5.2) 

Therefore let 5' G Fag(T z (H(u)), A) and 6% G Eig(T z (H(u)), A) x such that 5 = 5' + 5%. Then 
we obtain 

(E Tz (h(oj))({^})So,So) = (E Tz{H{uj)) ({\})5' ,5o) + (E Tz{h{lu)) ({X})6q,So) = (5' ,5 ) 
and with the above equivalence we get 

(6' ,6o) = (T_ z (5' ),T_ z (6 )) 

= (E HH ({\})T_ Z (5' ),T„ Z (5 )) + (E H(w) ({X})T. z (5^,T. z (5 )) 
= (E H(w) ({X})5 z ,5 z ), 



which implies ([52]). Applying (gTQ, Q and the fact T x (H(oj)) = H(T x (u)) leads to 

E <^fl(r.( w ))({A})^o,5o>. 



1t(XA n ^H( w )({A})) 1 



Finally we use Lemma 13.31 and Theorem l5.1l to obtain the existence of a set SlCfiof measure 
one such that for each to € ft we have 

lim r:— ; = / {E H ^){{X})So,S )dP{u) 

n ^°° |A n | Jq 

which was to prove. ■ 

The following fact is taken from |LV09j 

Lemma 5.3. Let r > 0, A C 7L d and U C £ 2 (A) be given and denote by U r the subspace of U 
consisting of all functions which vanish on d r (A). Then 

< dim([/) - dim(£/ r ) < |d r (A)|. 

Proof. Let P : U — >• £ 2 (<9 r (A)) be the natural projection with (Pip)(x) = tp(x) for all x G 
d r (A). Then we have 

< dim(J7) - dim(kerP) = dim(ranP) < \d r (A)\, 
which proves the claim as ker P = U r . ■ 

For given uj G 0, R G N and Q <Z Z d finite, let L^(R,Q) denote the number of e G -E 
with 6 e (cj) = 1 which are of length not less than R and incident to some vertex in Q, i.e. 

L^\R,Q) := \{{x,y} G E \ b {x>y} (oj) = 1, d(x, y) > R and {*, y} n Q ^ 0} | . (5.3) 

Let (A n ) be the sequence of cubes given as in (|3.12p . We chose a function R : N — > N such 
that 

\d R ^A n \ 

lim R(n) = oo and lim — — — : — = (5-4) 



n—too n— too 
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and set 

4 W) :=L^\R(n),A n ). (5.5) 

Beside this we set for R > 

xeZ d ,d(0,x)>R 

and for n € No 

e (n) := e R(n) as well as 5(n) := (2n + l)" d/4 . (5.6) 
Note as p € £ 1 (Z d ) we have by the definition of R(n) that 

lim e(n) = lim 5(n) = 0. 

n— >-oo n— >oo 

The next result estimates the probability that the number of long edges is large. 

Lemma 5.4. Let (Q,A,¥) be given as above. Then the following holds: 

(a) There exist constants Ro € N and 5 > such that for all < 5 < 5, all R > Rq and all 
finite QCZ d 

S 2 \Q\' 



F(LM(R,Q) > \Q\(e R + S)) <exp 



(b) Let R-.n be as in ([S3]) and L^ = L^(R(n),A n ). Then there exists a set Cl C Q 
of full measure such that for each ui £ f2 there exists no(ui) with 

4° < \AnMn) + 6(n)) (n>n (u;)). 

Proof. The proof of part (a) is to be found in |Schl2j . It is basically an application of a 
Bernstein inequality. Let us prove part (b). Therefore consider the events 

4:= 4") > \A n \(e(n) + 8(n))} . 

Then part (a) shows that for n large enough we have 

¥(A n ) < exp (-,5(n) 2 |A n |/4) = exp(-(2n + l) d ' 2 /A) , 
which clearly gives XmeN^(^«) < 00 • By the Lemma of Borel Cantelli we have 

P (lim sup A n ) = 

which implies the claim of part (b). ■ 

We use Lemmas 15.21 15.31 and 15.41 to obtain a result similar as Lemma 6.2 in |LV09] 

Theorem 5.5. Let (Q,A, P) be the probability space, H the randomly weighted Hamiltonian, 
(A n ) be the sequence of cubes and pffl the measures associated to the eigenvalue counting 
functions given as before. Then there exists a set £1 C Q of full measure such that for all 
well and oil A £ 1 we have 

lim ^lMl = E{(E H ({\})5 , 5 )}. 

n->oo |A n | 
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Proof. Let C Q be a set of full measure such that the results of Lemma 15.21 and of Lemma 
EH (b) hold for all uj G Cl. We fix some cj € Cl and AgM. With the function R : N ->■ N given 
in (|5.4p we set 

# := G l\l d ) | (fi"(o;) - A)u = and suppt- C A n „« (n) } , := dimV^. 

Note that Vn consists of the elements i An v, where v G ^ 2 (A n ) satisfying v = on A n \ 

(p An H(uj)i An -\)v = Q and ^ (H x ^(u) - \5 x {y))v{y) = (5.7) 

for all x € with x ~ A n _#( n ) . 

We consider the following difference 

\p%>\{\}) - Tr( XAn E HH )\ < \p^({\}) - DM\ + |DM " Tr(xA„^ H )| (5.8) 

and treat the two summands on the right hand side separately. Let us estimate the first one. 
Consider therefore the sets 

Ujfi := {u G £ 2 (A n ) | (p An H(u)i An - X)u = 0} 

and 

U n,R = i U e Un I U = ° n An \ A n-R(n)} ■ 

Then clearly, p ( n ] ({\}) = dim(C/^ ) ) > dim(K M ) and 

dim(E$>) - dim(y^) < \{y G A c n | y Z A n „ R(n) }| < L^(R(n),A n ) = (5.9) 

where we used the definition (|5.3p . The application of Lemma 15.31 gives 

< p^({A}) - DM = dim(C/M) - dim(yH) < dim(^)) - dim(C^) + 4 w) 

< l^^l + LH. (5.10) 

Now we estimate the second summand in (|5.8p . Therefore let «i, . . . , u n (w) be an orthonormal 

basis (ONB) of K H and let Vi, i G / be an ONB of the orthogonal complement of Vn"^ in the 
space Eig(iT(cj), A). Furthermore let Vj, j G J be an ONB of Eig(iT(a;), A)^. Then we have 

Tr( XAn EU{\})) = (XA n ^({A})^,^) + ^(XA n ^({A})^,^) + ^(xA„^({A})^,^) 
i=i ie/ ieJ 

D (-) 

i=l 

which gives < Tr(x An ^({A})). Now let m, i e I be an ONB of 

:= ran( XAn ^({A})) 
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and uj, j G J be an ONB of (U^) 1 . Then, using Cauchy Schwarz inequality, we obtain 

{XAn E oj(W) u i> u i) < IIXA n ^({A})«i||||«i|| < 1 and (xA„E UJ ({\})u j: u j } = 
for all i G I and all j G J. This gives 
4 W) < MxaM{\})) = Yl {XKM{X})u hUi ) + (XAMiX})^,^) < dim(C7H). 

iG/ j€.J 

(5.11) 

where we used dim(J7 n ) = |7|. As before we denote by U^r * ne subset of £/n consisting of 
the elements which vanish outside of A„_r. Therefore we have 



{xA n ^ | « € ^ 2 (Z d ), -A)« = 0,^0on d R ^A n } . (5.12) 



In the next step we define a set ^ 5 by dropping conditions in (|5. 12|) , in the following 
way 



U { :1 :=<XA n v 



G £ 2 (Z d ), {H x , y {u) - X5 x (y))v(y) = for all x G Z^\v = on c^A r 



SXA n V 



G ^ 2 (Z d ), ^ (flx,»(w) - \5 x (y))v(y) = for all xeZ<?\v = on ^A n L 
y&A n J 



where 



4 w) = A n _ fl(n) \ {x G A n _ fi(n) | x ~ A^}. 

Here we used that for all x G Zn and y G A^ we have H xy (uS) = 0. 

Comparing this representation of with the representation in (j5.T|) . we realize 

that they differ in at most 2Ln + |(?^^ n ^A n | conditions. As each of these conditions may 
change the dimension at most by one, we get 

dim(^) < dim(^) < DM + 2LM + |^A n |. (5.13) 

Applying (|5. 1 1 1) . Lemma I5T31 and ()5. 13|) gives 

< Tr( XAn ^({A})) - DH < dim(t^)) - < dim(t^) - D^ + |^A n | 

< 2\d R ^A n \ (5.14) 

In the last step we apply Lemma 15.21 then we combine the estimates for the two summands 
given in (|5.10p and (|5.14|) and finally use part (b) of Lemma 15.41 to obtain 

r ^(W) w/ P AYi i- U^«A})-Tr( XAn ^ H )| 
hm — — E((D if(a; )({A})do > <)o>) = hm 

n— »oo 



111 



I A, 



< lim 

n— voo 



I An | 

3|^)A n |+3L M 



< 3 lim 

n— ¥oo 



I An | 

'|d R (")A n | 

I A n | 



+ e(n) + <5(n) = 0. 



Here we used the definitions of R(n), e(n) and 5(n) in (|5.4p and ()5.6 
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Remark 5.6. (a) Let us stress the fact that proof of Theorem 15.51 does not contain any 
probabilistic argument. We show the claimed convergence for any fixed choice of A G R 
and well, where is a set given rather explicitly by Lemmas 15.21 and 15.41 

(b) Furthermore the proof gives an explicit error-term on finite scales. To be precise we 
have for any n G N, A G R and oj G f2 

\pfr\{\}) - Tr( XAn E H(u)) )\ < 3\d R ^A n \ + 3LM 

where L ( n ] = Z>)(i2(n), A n ) as in (EDI . 

The following result is essentially standard and has been used in the present context 
already in |LV09] . It shows that weak convergence of measures plus convergence of the 
measures at each point implies uniform convergence. 

Lemma 5.7. Let p be a probability measure on R and let (p n ) be as sequence of bounded 
measures on R which weakly converge to p and fulfill 

hm p n ({\}) = p({\}) 

n— >oo 

for all A G R. Then the distribution functions F n : R — > R, := Pn((~ °°>M) converge 

with respect to supremum norm to the distribution function F : R — > R, -F(A) := p((—oo, A]). 

The proof of the main theorem, already stated in Section [2] is now basically a combination 
of the previous lemmas. 

Theorem 12.11 Let (Q,,A,F) be the probability space, H the randomly weighted Hamiltonian 
and (A n ) be the sequence of cubes given as before. Furthermore let F£ be the distribution 
function given by (f3TT3|) and let F : R ->• R, F(A) := E{(E H ((-oo, \])S ,8 )}. Then there 
exists a set HCJ! of full measure such that for all uj G Q we have 

tt^t ->■ F as n ->• oo 

l A n| 

with respect to supremum norm. 

Proof. Let p, p { n ] : B(R) -> [0, 1] be the measures associated to the distribution functions 
F respectively F£. Then obviously p is a probability measure and the measures p^ are 
bounded. As shown in Corollary 14.21 there exists a set fii C with P(fii) = 1 such that for 
all a; G Qi the measure p is the weak limit of p^h \ Furthermore we have by Theorem l5.5l a set 
^ 2 C O with P(fi 2 ) = 1 such that for all w G tt 2 and all A G R one has p^ ({A}) p({A}). 
Therefore Lemma 15.71 yields the uniform convergence of the distribution functions for all 

w s Hi n fi 2 - ■ 
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